Magnetogenesis in bouncing cosmology 
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We consider the process of magnetogenesis in the context of nonsingular bounce cosmology. We 
show that large primordial magnetic fields can be generated during contraction without encounter- 
ing strong coupling and backreaction issues. The fields may seed large-scale magnetic fields with 
observationally interesting strengths. This result leads to a theoretical constraint on the relation 
of the energy scale of the bounce cosmology to the number of effective e-folding of the contracting 
phase in the case of scale invariance for the power spectrum of primordial magnetic fields. We show 
that this constraint can be satisfied in a sizable region of the parameter space for the nonsingular 
bounce cosmology. 


PACS numbers: 98.80.Cq 


I. INTRODUCTION 


Astronomical observations indicate the existence of magnetic fields throughout the observable universe-from the 
length scales of stars up to the scales of galactic clusters. The strength of these fields in galaxies and clusters is 
on the order of u Gauss. Magnetic fields can be produced from the amplifications of initial seeds via some physical 
processes during structure formation (such as adiabatic compression and turbulent shock flows). Observations suggest 
the strengths of initial seed fields should be greater than 10713 Gauss on a comoving scale larger than 1 Mpc [I B]. 
If the seed magnetic fields originated in the very early universe, before the surface of last scattering, they may effect 
the temperature and polarization spectra of the cosmic microwave background (CMB) anisotropies since their energy- 
momentum tensor can source the corresponding scalar, vector, and tensor cosmological perturbations. Therefore, CMB 
observations can provide experimental constrains on the properties of primordial magnetic fields [BH9]. According to 
the latest Planck 2015 data, primordial magnetic field strengths are required to be less than a few 107° gauss at the 
length scale of 1 Mpc [0]. 

Discovering the underlying mechanism responsible for generating primordial magnetic seed fields is an active area 
of cosmological research. In the literature, the origin of primordial magnetic fields has been widely discussed, see 
Refs. [IIHI3] for relevant reviews. Amongst these mechanisms, one possibility is that the seed fields could be produced 
during cosmic phase transitions such as the electroweak or QCD phase transitions. This approach, however, suffers 
from a manifest issue that the correlation lengths are too small when compared with observations. Consequently, one 
may investigate the possibility of producing seed fields at an even earlier time, such as during the inflationary epoch, 
in which case the process is known as inflationary magnetogenesis. Since the theory of electrodynamics is conformally 
invariant, and the energy density of the background universe is nearly constant, inflationary magnetogenesis requires 
the breaking of conformal invariance by introducing couplings of the electromagnetic fields to other sectors such 
as gravity [I4HI8], inflaton fields [19}{24|, pseudo scalar fields 25H27], or other fields [28430]. However, these 
models typically produce seed fields that have too small strengths or suffer from other conceptual issues, such as the 
backreaction problem that the energy density of the generated electromagnetic fields spoil the inflationary background 
dynamics, the strong coupling problem that the coupling of electromagnetism (EM) may become too large in the 
primordial stage. Additionally, there exists the curvature perturbation problem wherein the curvature perturbations 
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seeded by the generated electromagnetic fields are too large to be in agreement with the CMB data 36]. Recently, 
two of the current authors proposed a generalized model within the inflationary scenario in which the above issues 
were avoided under a certain fine tuning [37]. 

One possible alternative to the inflationary paradigm, is that our observed universe may have evolved from a pre- 
viously contracting phase and then experienced a nonsingular bouncing phase at very early times. This is known as a 
bouncing scenario [38H43]. In particular, it has been shown that dynamical bounce models under certain contracting 
phases, such as the matter contraction [44H46] or the ekpyrotic phase [42] [47] [48], a nearly scale invariant power 
spectrum of primordial fluctuations could be achieved in order to explain the present cosmological observations. Var- 
ious model constructions and related perturbation analyses were extensively performed in the literature, namely, the 
quintom bounce [49H51], the Lee-Wick bounce [52], the Horava-Lifshitz gravity bounce [53H55], the f(T) teleparallel 
bounce [56H58], the ghost condensate bounce [59], the Galileon bounce [60] [61], the matter-ekpyrotic bounce [62H64], 
the fermionic bounce [65] [66], etc. (see, e.g. [44] [67] for recent reviews). In general, it was shown that on length 
scales larger than the time scale of the nonsingular bouncing phase, primordial cosmological perturbations remain 
almost unchanged through the bounce [69]. The bouncing phase embedded in inflationary scenario was recently 
considered to explain both the large hemispherical power asymmetry and the power deficit on large angular scales in 
cosmic microwave background radiation [70] [71] (see [72] [73] for the bouncing inflation). 

In this paper we investigate the process of magnetogenesis within the frame of the nonsingular bounce cosmology. 
The magentogenesis in bouncing universe was once investigated in where a specific form of the gauge field 
coupling was considered. Different from the case in inflationary scenario, the bounce cosmology suggests that our 
universe had a sufficiently long period of contraction and then experience a smooth bouncing phase that connected the 
contracting phase with the observed expanding one. Correspondingly, the energy density of the background universe 
cannot be approximated to be a constant as in the inflationary scenario. It is notoriously difficult to construct bouncing 
models in General Relativity in a flat universe due to the required violation of the null energy condition (NEC). In 
general, NEC violation in bouncing cosmology is accompanied by a number of stability problems. If the matter 
Lagrangian is simply a general function of a scalar field and its derivative (as in k-essence) one cannot transition into 
an NEC violating phase and hence a bounce is not possible [74] [75]. A construction of a bounce which was manifestly 
free from gradient and ghost instabilities before thermal expansion and which occurs on trajectories of a non-vanishing 
measure was achieved in [61]. As a specific model for bouncing cosmology we choose [62] to study magnetogenesis; 
however, our findings are applicable to general bouncing models. In this paper we show that a bouncing universe 
can realize magnetogenesis in an efficient way due to the following three conditions. Firstly, as previously mentioned, 
there exists the large scale problem for the generation of primordial magnetic field. This issue is nicely addressed 
within the bouncing scenario since in the contracting phase the wavelength of the seed fields exit the Hubble horizon. 
Secondly, the scale factor in the contraction phase is decreasing and the power spectrum of the magnetic field becomes 
dense instead of being diluted as in inflation. Thirdly, the energy density of the background universe also becomes 
dense in the contracting phase, and therefore the backreaction problem can be avoided within a reasonable tuning of 
model parameters. 

The paper is organized as follows. In Section II we introduce the basic formulae of the model to realize magneto- 
genesis in bouncing cosmology. In Section III we briefly review the scenario of bounce cosmology by considering a 
specific model but keep the background evolution general. We then analyze the generation of primordial seed fields 
in Section IV and discuss constraints from the strong coupling issue. The backreaction issue of primordial magnetic 
fields is investigated in Section V within the bouncing scenario. We summarize our results in Section VI. 


I. THE MODEL OF BOUNCE MAGNETOGENESIS 


To begin, we consider the following action for scalar field ¢ and electromagnetic field A 
Einstein gravity 


u, minimally coupled to 


S= ftev ER + Loon = POF , (1) 


where Lpounce is the Lagrangian density responsible for the nonsingular bouncing background which will be introduced 
in detail in Section III, and the field strength F,,,, is defined in terms of the U(1) gauge field in the familiar way 
F v = ô Ay — 0,A,. Unless explicitly stated, we will work with natural units in which c = h = Mp = 1. 

Moreover, the standard Lagrangian of the electromagnetic field can be recovered when the coefficient f is taken to 
be unity. In order to achieve the amplification of the magnetic fields, one expects the coefficient f to be a function 
of the background scalar field ¢ and hence can evolve along with the cosmic time so that the conformal symmetry is 
violated under this time evolution. 


Varying the above action with respect to the electromagnetic field yields the corresponding equation of motion 
as follows, 


V, [ -F FY] =0. (2) 
We consider an observationally favored, spatially flat Friedmann-Robertson-Walker metric, 
ds? = a*(n)(—dn? + dx”), (3) 


where a(n) is the scale factor and 7 is the conformal time. 
We choose the Coulomb gauge: Ag = 0, 0;A’ = 0, and expand A; in the Fourier space as follows, 


3 
A;(n, x) = 5 Eo (ci o(k)ak,oAr(n)e™* + h.c.) ; (4) 


g=1,2 


where €;,,(k) are two orthonormal and transverse polarization vectors. The operators ak, and al > are the annihi- 
lation and creation operators that satisfy the standard commutation relations |ak,c, al, ot] = 900'0(k — k’) as well as 
[@k,05 Uk’,o’] = lai o al cl =0. 

Within the metric (3), the Fourier mode A; with a fixed conformal wave number k satisfies the following equation 
of motion, 


Al + °F A + kA, =0, (5) 
with the normalization condition 
An MAG O) = Au (MAL) = =z (6) 
being satisfied. In terms of a new variable vg = f Ax, the equation of motion can be reformulated as: 
whe- £) vk = 0. (7) 


From this equation, one can explicitly read that the evolution of the magnetic field, which is depicted by the gauge 
field A;, depends on the form of the coupling coefficient f and the background evolution. In the following section we 
investigate the background of a bouncing universe. 


II. BOUNCE COSMOLOGY 


Nonsingular bounce cosmologies appear in many theoretical frameworks where the gravitational sector is modified, 
by making use of Null Energy Condition violating matter fields, or in a universe with non-flat spatial geometry (see 
e.g. (77). However, it is important to note that the BKL instability could appear in the contracting phase since 
the back-reaction of anisotropies would dominate over the regular dust and radiation densities, unless one finely tunes 
the initial conditions to be nearly perfectly isotropic. This issue can be avoided if a period of the ekpyrotic phase 
is introduced before the bounce [78] where the unwanted anisotropies are depressed in the presence of an ekpyrotic 
scalar field. 

Recently, it was shown that nonsingular bouncing cosmology with regular matter contraction can be achieved with 
an era of ekpyrotic contraction by introducing a scalar field with a Horndeski-type, non-standard kinetic term and 
a negative exponential potential [62]. One could include a regular dust component and then assume the universe 
began in a state of matter-dominated contraction thus combining the matter bounce with the ekpyrotic scenario. 
The backreaction of anisotropies is manageable throughout the entire cosmological evolution of the matter-ekpyrotic 
bounce model, including at the bounce point [79], and hence, the BKL instability that arises for a large family of non- 
singular bounce models, as pointed out in BI], is nicely resolved. A concrete realization of the matter-ekpyrotic 
bounce was constructed in that involves two matter fields. This effective field theory model of a non-singular 
bounce can be embedded into a supersymmetric version [82], or into loop quantum cosmology [64] (also see [83H85] 
for different cosmological scenarios). 


In the present study, we follow the original bounce model within the context of effective field approach that was 
proposed in and its cosmological implications to magnetogenesis. The lagrangian can be written as 


Lyounce = —K(¢, X) a G(¢ġ, X) Q, (8) 


where K and G are functions of the dimensionless scalar field ¢ and its canonical kinetic term X = 0,¢0"¢/2. The 
d’Alembertian operator is given by Dd = g#”V,V_¢. 
Following [62] the function K is chosen as 


K(¢, X) = (1 — g(¢))X + BX? - V(¢), (9) 


and G as 
G(X) = 7X, (10) 
where 3, y are the related model parameters. For a phase of ekpyrotic contraction, the form of the potential V (¢) is 


taken to be 
2Vo 


where Vo, q ,by are model parameters as well. In order to obtain a violation of the Null Energy Condition near the 
end of the contraction phase, we take g(@) to be 


290 
gd) = (12) 
e VE 4 cdot 
where go,bg and p are model parameters whose reasonable range can be found in reference [62]. 
The evolution equation of the background dynamics are evaluated explicitly in this model, and are well studied in 
the literature. We do not repeat this study here, but rather provide a summary of the three stages of evolution most 
relevant to the process of magnetogenesis: 


Vig) = (11) 


e the ekpyrotic contraction, 
e the bouncing phase, 


e the fast-roll phase. 


Note that, compared to the first phase, the effects of the latter two are rather limited based on the analyses in [68], 
and hence, we treat them as instantaneous for simplicity. Moreover, we are interested in the contraction phase, during 
which the wavelength of the fields exits the Hubble horizon since this length scale contracts much faster than the 
scale factor. It is important to notice that during the same stage the energy density of the scalar field ¢ evolves as a 
perfect fluid with an almost constant equation of state w. 

During this phase ¢ < —1 so that g — 0 and ¢ < 1. The Lagrangian for ¢ approaches the standard canonical 


form. There is an attractor solution for œ which yields an effective equation of state, 
1+ 2 (13) 
wre —. 
3q 


If q has a critical value, the equation of state w can be sufficiently large and the energy density of ọ will dominate 
over matter and radiation during contraction. The background in this case can be described as 


a a7 
a x (fje —)™4, (14) 
where the subscript e stands for the end of the contraction phase, 7). = Ne — UDa T Ne is the end of the contraction 


phase and 7 = 0 is the midpoint of the bounce scenario. The above description is valid until we reach ne. For 
simplicity, we define g = q/(1 — q) and the scale factor can be expressed as 


~ \—-G 
ol) =a (=r) emt (15) 
where ae is the scale factor at the end of the contraction phase. 

After the moment ne, the universe can evolve through the nonsingular bouncing phase smoothly. During this 
phase, the evolution of the Hubble parameter H can be approximated as a linear function of the cosmic time t 
[63]. Accordingly, the scale factor behaves as a ~ exp(t?) roughly. In this case, one can read that a nonsingular 
bounce can occur when the scale factor reaches the minimal value. 

In the following section we perform the analysis of the generation of the primordial magnetic fields in the afore- 
mentioned cosmological scenario. 


IV. BOUNCE MAGNETOGENESIS 


After the bounce we wish to recover classical electromagnetism in order to connected with the observed big bang 
history smoothly, and therefore we require f — 1 along the background evolution. If in the contraction phase f is 
much smaller than unity we will face the so-called the strong coupling problem. A natural way to solve such a problem 
is to require f? > 1 during the contraction phase and it becomes unity at the end moment of the contracting phase 
Ne, Which is also the beginning of the bouncing phase. This can be done by assuming that, during the contracting 
phase when 7 < Ne, f is a power-law function of the scale factor via 


j=} (+). (16) 


Qe 


where fe ~ O(1). Since in the contraction phase the scale factor decreases, we take n > 0 so that f > 1 is fulfilled. 
During the contraction phase, the Hubble horizon decreases faster than the scale factor so that the modes with larger 
wavelengths will exited the Hubble horizon earlier. 

For short wavelength modes, or k? > f”/f, the solution to Eq. (7) is approximated as 


= —tky 
v = ——e ; 17 
while long wavelength modes (k? < f”/f) have the general solution, 
dn 
p(n) = Cif + Caf P (18) 


where C and C2 are integration constants that can be fixed by matching this solution to (17). With the background 
a (ñe — n)! and integrating the second term we find 


ug (7) = Cia” + Cie ta, (19) 


If —1 — n + 1/q < n, the second term will dominate over the first in the contraction phase. Using the conjunction 
condition 


1 , 
Cra, Se | a (20) 


where a, stands for the scale factor when the corresponding wavelength exits the Hubble horizon a| H| = k (which 
is also equivalent to k(fje — nk) = q) one derives 


~ 1 
l on+i-2 ae He) E00) 
C2 = Ja” j (- k ) : (21) 
We arrive at the complete form for vk, 
1 He q(n+1—-+) —n-1+4 
v = — | -2 = (22) 
V2k k de 


Actually, with the scale factor given by Eq. and the coupling function given by Eq. (16), the analytical solution 
to Eq. is a linear combination of Bessel functions [22]. The two integration constants can be fixed by imposing 
the Bunch-Davies Vacuum condition as (—kn) — oo. For modes well outside the Hubble radius (k? < f/f), 
one can get the same solution as (22) in the case of —1 — n + 1/q < n. Therefore, the gauge field is found to be 


1 aH ã(n+1— 4) ä —2n-1+4 
A, = vz, /f = ——— | -< — . 23 
mad Al k ) 3 a 
Correspondingly, the power spectrum of the primordial magnetic fields is evaluated as 
2q(n+3)—6 An—642 
k| Ak (n)|? 1 aeHe Tq a nd 2 5 
k.n) = = as H 24 
Pal 1) 5 2r2a4 Qn? f2 k de e? ( ) 


o=1,2 


with the spectral index ng = (6 — 2(3 + n)q)/(1 — q). Note that the spectrum is scale-invariant for n = —3 + 3/q. At 
the end of the contraction phase we have 
2q(n+3)—6 
=a 


1 aH, 
Pak, ae) = Qn? f2 (- k ) He (25) 


As argued previously, the effects of the (fast) bounce are limited and hence we treat this phase instantaneously. 
Note that the vector field has become the standard electromagnetic field right before the bounce at the moment ne. 
After that, the magnetic fields decay with the square of scale factor. The strength of the magnetic field observed 
today is given by: 


a(n+3)—3 
1 aeH, Iza Q? 
By =(P k, 1/2 = _ Mette £ H2, 26 
0 ( B( ao)) 2r fe k az e ( ) 
where ag is the scale factor today. In the scale-invariant case we have 
1 a2 
Bo = e 27 
g V2 fe ° ag ( ) 


We must determine the amount the scale factor changes from the end of the contraction until today. Using the 
conservation of entropy we find 


na (28) 


We treat the energy density of the electric field as the energy density of radiation at the end of the contraction 
phase because it dominates over that of magnetic fields. The power spectrum of the electric fields takes the form of 


P k —_ k?| Aj, (n)|? a (2n ie 1 = a aeHe ag) a “ee Ht (29) 
B( ,1) = 2. Ir2a4 Qn? f2 k de et 
The ratio of Pg to Pg can be written as 
PE 2/ 4A a= aeHe : 
_“=(2 ae ms 30 
nti- (E) "(3% ), (30) 


from which we can see that Pz > Pe for the long wavelengthes beyond the Hubble horizon. Therefore, we can 
compute the energy density of the electric part as the total energy density of the electromagnetic fields 


eH, 

etle dk 

PEM = rI PET 
aiH; 


(—1+ q+ 2ngq)?(1 — q) (ae He) ICHA / a N7 aH a “S82 
Qelle =g — (illi =g , 
27?q?(4 — 2(2 + n)q) 


(31) 


Qe 


where a; is the scale factor at the moment when the largest scale observed today first exited the Hubble horizon. 
Since at the end of the contraction phase a;H; < aeHe in the case of the scale-invariant power spectrum, the first 
term on the right-hand side of equations can be ignored. Thus, we have 


25(1 — 9)*aeHe 


—2 
272q2 (a; H;) ’ (32) 


PEM(Ge) ~ 


From equations and (28), we get 


ge E vV2Tênq aiHi 


H . 33 
Cag =5(1—q) aeHe (33) 
The strength of the magnetic fields today is then given by 
Te iA; Te 
Bo o1 2t exp(—N), (34) 


T 5fe(1 — q) acHe ~ 5fe(1 — q) 


where N is the effective e-folding number of the contraction phase. Following Ref. [79], it is limited by 


aie) (=) (35) 


where Qpert is the amplitude of the observed anisotropy in CMB. This leads to 


Tod o 
Bo < 5f =q) pert ’ (36) 


Choosing fe = 1 and q = 0.1, and using Tp = 1.9 x 1078? in Planck units and Qyert = 1075, we have the upper bound 
of the strength of the magnetic fields today: 


Bo < 1071! Gauss. (37) 


Note that, we need to ensure that the contracting phase proceeds for a sufficiently long period as discussed above, 
and the condition that —1 — n + 1/q < n has to be satisfied in the scale-invariant case. The combination of these two 
conditions then yields the upper bound of Bo as a function of the model parameter q. In Fig. |1| we demonstrate the 
relation between the upper bound of Bo and the parameter q for specific parameter choices. 


log,.(Bo) 
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FIG. 1: Dependence of the upper bound of the magnetic fields Bo on the model parameter q in the bounce model under 
consideration. Other model parameters have been taken to be the same as those provided in Ref. [79]. 


V. CONSTRAINTS FROM THE BACKREACTION 


In the above section we performed the analyses of the magnetogenesis within bounce cosmology and showed that 
the corresponding results can satisfy the observational constraints which gives rise to an upper bound for the produced 
magnetic fields. In this section we further check the conditions that the energy density of the electromagnetic fields 
does not spoil the evolution of the background dynamics, which requires 


PEM < Pọ, (38) 


where pọ stands for the energy density of the background’. Since ¢ dominates over other components at late times, 
we may compare it with the contribution of the electromagnetic fields directly for simplicity. From we find 


pem xa 2"-4 x az a, (39) 
in the case of scale invariance. Moreover, the energy density of the background scaler field evolves as 


Po x GT yami, (40) 


1 Note that, if there exist degrees of freedom other than the background field ¢, such as the radiation, a second matter field, or even the 
cosmic anisotropy, the stability of the bouncing phase can be well under control in certain parameter space of the bounce model, as has 
been comprehensively studied in . From [79], one can see that a nonsingular bounce can be obtained when the total energy density is 
vanishing. In this case, the background dynamics of the universe is dominated by the pressure of the field ¢ instead of energy densities. 
As a result, the requirement pgm < pg remains valid near the bounce. 
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The above two equations imply that the energy density of EM grow even faster than that of ¢. We require that at 
the end of the contraction phase ¢ dominates, that is to say, 


pem(ae) < H2M3,, (41) 


where we recover the Planck energy scale unit Mp, to balance the units on both sides of the above equation. This 


gives 


25(1 — q)?a2H? 
nA HNM <b (42) 
27?q? (a; Hi) Mi 


which leads to 


|He| _ vV2rqexp(-N) 


My ~< 50-9) 12 


As a result, the inequality yields a theoretical constraint on the relation of the effective e-foldings N to the energy 
scale of ¢ at the end of contraction phase. 

In order to have a clearer picture of this theoretical constraint, we numerically depict the parameter space that is 
allowed for a well behaved bounce model in Fig. |2} We find this theoretical constraint is easily satisfied in a sizable 
region of the parameter space of nonsingular bounce cosmologies. 
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FIG. 2: The theoretical constraint of the energy scale of the bounce cosmology |He|/M,: along with the model parameter q. 
The green shadow area satisfies the theoretical constraint and hence bounce models within this regime are well behaved. 
We take the effective e-folding number to be N = 10 as an example. 


VI. CONCLUSION 


In this paper we showed it is possible to produce appreciable primordial magnetic fields within the context of 
nonsingular bounce cosmology. We adopted the specific model developed in [62] but the results derived in our study 
are applicable to more general bouncing models as the main analyses only depend on the phase of ekpyrotic contraction. 
In this model, we found that primordial magnetic field with potentially sizable strength could be produced without 
troublesome strong coupling and backreaction effects, which plague the process of magnetogenesis in other very early 
universe paradigms such as inflation. By focusing our interest on the nearly scale-invariance power spectrum of 
primordial magnetic fields, we derived an explicit observational constraint upon the relation between the magnetic 
field that may be observed today and the model parameter that is associated with the effective equation of state during 
the contracting phase. Moreover, by examining the limit of the backreaction, we obtained a theoretical constraint 
on the relation of the energy scale of the bounce cosmology with respect to the effective e-folding numbers of the 
contracting phase in the case of scale invariance. We showed that this theoretical constraint can be satisfied in a 
sizable parameter space of bounce cosmologies. 
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